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We present new accurate measurements of the differential cross section a (8) and the proton 
analyzing power A y for proton- 3 He elastic scattering at various energies. A supersonic gas jet target 
has been employed to obtain these low energy cross section measurements. The cr{6) distributions 
have been measured at E p = 0.99, 1.59, 2.24, 3.11, and 4.02 MeV. Full angular distributions of 
A y have been measured at E v = 1.60, 2.25, 3.13, and 4.05 MeV. This set of high-precision data is 
compared to four-body variational calculations employing realistic nucleon-nucleon (NN) and three- 
nucleon (37V) interactions. For the unpolarized cross section the agreement between the theoretical 
calculation and data is good when a 37V potential is used. The comparison between the calculated 
and measured proton analyzing powers reveals discrepancies of approximately 50% at the maximum 
of each distribution. This is analogous to the existing " A y Puzzle" known for the past 20 years in 
nucleon-deuteron elastic scattering. 

PACS numbers: 21.45.+V, 21.30.-x, 24.70.+S, 25.40.Cm 



I. INTRODUCTION 



The study of the four nucleon (47V) system is interesting from a number of different perspectives. First of all, many 
reactions involving four nucleons, like 2 H(d,p) 3 H, 2 H(d,n) 3 He, or p + 3 He — > 4 He + v c + c + (the hep process), are of 
extreme astrophysical interest, as they play important roles in solar models or in big-bang nucleosynthesis (BBN). 
The hep process, for instance, is the source of the highest energy neutrinos from the Sun. Moreover, 47V systems have 
become increasingly important as testing grounds for models of the nuclear force. While not the nuclear-structure 
"imbroglio" of heavy nuclei, the A = 4 system is the simplest system that presents the complexity — thresholds and 
resonances — that characterize nuclear systems, and therefore is a very good testing ground of modern few-body 
techniques Q. Similarly, since the 47V bound state is (to a very good approximation) a ( J 17 , T) — (0 + , 0) state, it is a 
good "laboratory" for the study of the strange quark components of the nucleon via parity- violating electron-scattering 
experiments 0- 

The theoretical description of A = 4 systems still constitutes a challenging problem from the standpoint of nuclear 
few-body theory. Only recently, with the near-constant increase in computing power and the development of new 
numerical methods, has the study of the a-particle bound state reached a satisfactory level of accuracy; the 47V 
bound-state has been calculated to a few tenths of keV The study of 47V scattering states, on the other hand, 
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is less satisfactorily developed. The same increases in computational power, however, have opened the possibility for 
accurate calculations of the 47V observables using realistic models for nucleon-nucleon (NN) and three-nucleon (37V) 
forces. These calculations have beenperformed mainly by means of the Faddeev-Yakubovsky (FY) approach 0,0,13 
and the Kohn variational principle |9l llflllll lL2lll3 | . In spite of this progress, some disagreements still exist between 
theoretical groups, as for example, in the calculation of the n- 3 H total cross section in the peak region (around 
Ben.. = 3 MeV) 0. 

To make matters worse, even in cases where the theoretical calculations agree, they are often strongly at variance 
with the experimental data. For example, the proton analyzing power A y in p- 3 He elastic scattering is underestimated 
by theory at the peak of the angular distribution by about 40% [10(. Other disagreements are discussed in Refs. 

Blllliip. 

The existing data in the literature for AN scattering are of lesser quality, when compared with the excellent and 
abundant data that exists for the NN and 37V (N-d scattering) systems. The intensely studied 4 He system is 
(unfortunately) very difficult to describe theoretically due to the presence of the bound state and many higher-energy 
resonances 16]. This leads us to investigate p- 3 He scattering states. 



A. n- 3 H Elastic Scattering 

This paper will focus on p- 3 He elastic scattering at low energies, which is simpler than the He system to investigate 
theoretically. However, the situation for the closely related n- 3 H elastic scattering case is worth briefly discussing 
first. The quantities of interest in n- 3 H zero-energy scattering are the zero-energy total cross section ctt(0) and the 
coherent scattering length a c . Experimentally, only the total cross section or has been measured with high precision 
for a large range of energies. The extrapolation of the measured ctt to zero energy is straightforward and the value 
obtained is ctt(O) = 1.70 ± 0.03 b ^?|]- The coherent scattering length has been measured by neutron- interferometry 
techniques. The most recent value reported in the literature is a c — 3.59 ± 0.02 fm 18]. An additional estimation 
of a c — 3.607 ± 0.017 fm has been obtained from p- 3 He data by using an approximate Coulomb-corrected R-matrix 
theory [l^. These values should be compared with those obtained theoretically: or(0) = 1.73 b and a c = 3.71 fm 9]. 

As already mentioned, at higher energies (specifically in the "peak" region) there exist sizable discrepancies between 
different theoretical predictions. The ctt calculations by Fonseca Q and Pfitzinger, Hofmann, and Hale 0] are i n 
good agreement with the experimental data, while the calculations by the Grenoble 0, ITil 12(1 ] and Pisa groups 0] 
are well below the data. The origin of this disagreement is still not clear. 



B. p- 3 Hc Elastic Scattering 

Let us consider now the situation for p- 3 He elastic scattering. The zero-energy quantities for this case are more 
difficult to evaluate. Approximate values of the triplet and singlet scattering lengths have been determined from 
effective range extrapolations [2l| to zero energy of data taken mostly above 1 MeV, and therefore suffer large 
uncertainties [ljj. This problem has been reconsidered recently by George and Knutson [l5|: this new phase-shift 
analysis (PSA) gave two possible sets of scattering-length values, both of which are at variance with the theoretical 
estimates (see the discussion in Ref. [la].) 

The world database of existing p- 3 He scattering data can be divided in three energy regions. There is a set of cross 
section and proton analyzing power A y measurements at very low energies, from E p — 0.3 to 1.0 MeV 22] . Another 
energy region is for E p = 1.0 to 4.0 MeV which includes the recent A y measurements at E p = 1.60 and 2.25 MeV 
|10|. However, these A y data are not very precise. The cross-section data in this region are similarly imprecise and 
sparse in number. They are also very old, being the very first published cross-section data for p- 3 He elastic scattering 
found in the literature j23|. The third group of measurements (for E p > 4.0 MeV) includes differential cross sections, 
proton and 3 He analyzing power measurements, and spin correlation coefficients of good precision. References for all 
these measurements can be found in Ref. pl|. 

The calculations performed so far for p- 3 He scattering have shown a glaring discrepancy between theory and 
experiment in the proton analyzing power [Tol Il3| . This discrepancy is very similar to the well known u A y Puzzle" 
in N-d scattering. This is a fairly o ld p roblem, already reported almost 30 years ago 0, |2^ in the case of n-d and 
later confirmed also in the p -d case |2(| |23 ■ All N-d theoretical calculations based on realistic NN potentials (even 
including 37V forces) underestimate the measured nucleon vector analy zing po wer A y by about 20 — 30 %. The same 
problem also occurs for the vector analyzi ng p ower of the deuteron iTn |26tl27| . while the tensor analyzing powers T20, 
T21 and T22 are reasonably well described |27ll28T]. The inclusion of standard models of the 37V force have little effect on 
calculations of these observables. To solve this puzzle, speculations about the deficiency of the NN potentials in 3 Pj 
waves (where the spectroscopic notation 2S+1 Lj has been adopted) have been suggested. Looking at NN scattering 
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data only, this possibility does not seem to be ruled out |29|,|30|. However, after taking theoretical constraints into 
account |3l|. this solution was considered unlikely, and has been confirmed by recent calculations which show that 
this puzzle is not solved even when new TV TV potentials derived from effective field theory are used 0, ■ 

Consequently, attention has focused on exotic 37V force terms not contemplated so far. The best solution seems to 
be given by the inclusion of a spin-orbit 37V force |3^| ; such a force has a negligible effect on the observables already 
well reproduced by standard 3TV potential models (such as the binding energies, N-d unpolarized cross sections and 
tensor analyzing powers) but it is very effective for solving — or at least for noticeably reducing — the discrepancy. 
Other explanations have been proposed in Refs. |3> HH HH US It should be noted that the current understanding 
of the 3N interaction is rather poor [38l l3p. and new terms in the 3TV interactions derived from chiral perturbation 
theory have very recently been proposed |4Ct l4l| . These new models have to be tested primarily in the 37V system, 
but the 47V system could also play an important role. In fact, 37V force effects are expected to be sizeable in the 47V 
system m f ne N-d scattering system 37V force effects are small [28| and masked in part by the contribution of 
the Coulomb potential 0,0]- Moreover, the N-d system is essentially a pure isospin T = 1/2 state. Tests of the 
T = 3/2 channel in any 37V force can only be satisfactorily performed in a 47V system such as p- 3 He. 

Seeking to explore potential three-nucleon force effects in a different system, and to investigate this new A = 4 
A y Puzzle, a series of proton- 3 He elastic scattering measurements have been made. Angular distributions of the 
differential cross section cr(9) and proton analyzing power A y have been measured at several energies below 5 MeV; 
the analyzing power experiments were performed with a gas-cell target and the cr(9) measurements were performed 
using a supersonic gas-jet target. The nominal proton energy E p at which each experiment was performed is listed 
in Table [HJ These energies were chosen to maximize the improvement in the low energy database, and because 
theoretical calculations are tractable in this region. 

In addition, we present new theoretical calculations for the p- 3 He system. They are performed via an expansion of 
the wave functions of the scattering states in terms of a hyperspherical harmonic (HH) basis and using the complex 
form of the Kohn variational principle |45l . These new calculations reach a much higher degree of accuracy than 
those performed previously using the correlated hyperspherical harmonic (CHH) functions fOj. In this paper we 
present calculations based on the Argonne v\% (AVI8) J47j 7V7V potential, which represents the 7V7V interaction in 
its full richness, with short-range repulsion, tensor and other non-central components and charge symmetry breaking 
terms. The calculations are performed without and with the inclusion of the Urbana IX (UIX) 37V force |48| . 

The remainder of this paper is organized as follows: in Section^ a brief description of the HH technique is reported; 
in Section lTTTl the experimental setup and methods are discussed; the cross section a(9) and proton analyzing power A y 
measurements are compared with the theoretical results in Section llVl finally, Section [V] is devoted to the conclusions 
and broader context of the present work. 



II. THE HH TECHNIQUE FOR SCATTERING STATES 

The wave function 9^^ Jz7r describing a p- 3 He scattering state with incoming orbital angular momentum L and 
channel spin S (S = 0, 1), total angular momentum J, and parity it = (— ) L can be written as 

\£jLS J J Z 7T ^LS J J Z 7T _|_ ^jLS J J Z 7T ^-j^ 

where \$/^ SJJ *' K vanishes in the limit of large intercluster separations, and hence describes the system in the region 
where the particles are close to each other and their mutual interactions are strong. On the other hand, \f/LS JJ ' 7T 
describes the relative motion of the two clusters in the asymptotic regions, where the p- 3 He interaction is negligible. 
In the asymptotic region the wave function ^SJJ,* 

reduces to JJz7T , which must therefore be the appropriate 
asymptotic solution of the Schrodinger equation. -qj 1 ^ SJJ - 7r ca n be decomposed as a linear combination of the following 
functions 

o± 1 V^Tr o^rini c - <\ (ti \ G L( r )c,qy l ) ,.F L (r) C ,qyi)\ 

where is the distance vector between the proton (particle i) and 3 He (particles jkl), q is the magnitude of the 
relative momentum between the two clusters, Si the spin state of particle i, and <p^ is the 3 He wave function. The 
total kinetic energy T c m . in the center of mass (cm.) system and the proton kinetic energy E p in the laboratory 
system are 

q 2 4 

T c .m. = ' g-^c.m. i (3) 
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where fj, = (3/4)M/v is the reduced mass (Mm is the nucleon mass.) Moreover, Fl and Gl are the regular and 
irregular Coulomb functions, respectively, with rj C — 2p,e 2 /q. The function fhiVi) = [1 — ex p(~72/i)] 2i+1 m @ nas 
been introduced to regularize Gl at small t/,, and fL(Vi) — > 1 as ?/i becomes large, thus not affecting the asymptotic 
behavior of ^f?^"^. Note that for large values of qyi, 



fL(yi)G L (r) C , qyi) ± iF L (r]c, qyi) -> exp ±i(gy, - Ltt/2 - 77c m(2gi/i) + ct l ) 



(4) 



and therefore, ^^SJJ (^lsjj ) describes in the asymptotic regions an outgoing (incoming) p- 3 He relative motion. 
Finally, 



LSJJ.v 



L'S' 



>LL 



s ss ,n 



LSJJ Z 



(5) 



where the parameters S^s L'S'io) are the S'-matrix elements which determine phase-shifts and (for coupled channels) 
mixing angles at the energy T c m .. Of course, the sum over L' and S' is over all values compatible with a given J and 

IT. 



parity. In particular, the sum over V is limited to include either even or odd values such that (— 1) 

c 



L' 



The "core" wave function JJ * 7T is expanded using the HH basis. For four equal mass particles, a suitable choice 



of the Jacobi vectors is 



tip 



x 2p 



(6) 



^-3p — *j M 



where p specifies a given permutation corresponding to the order i, j, k and m of the particles. By definition, the 
permutation p = 1 is chosen to correspond to the order 1, 2, 3 and 4. For a given choice of the Jacobi vectors, the 
hyperspherical coordinates are given by the so-called hyperradius p, defined by 



/> x /xl p + x\ p + xl p , (independent of p) , (7) 

and by a set of angular variables which in the Zernike and Brinkman |4flL representation are i) the polar angles 
Xi p = [0ip,4>i p ) of each Jacobi vector, and ii) the two additional "hyperspherical" angles 4>2p and 03 P defined as 

xip , *^3p 



COS< 



2p 



cos< 



'3p 



'lp 



l 2p 



(8) 



b 2p 



u 3p 



hereafter as tt p . A generic HH function is written as 



where the Xi P are the moduli of the Jacobi vector Xi. The set of angular variables x\ p , X2 P , Xz p , 4>2 P , and 4>z p is denoted 

(9) 



y 



K,A,M 

1,(2, (3, L2, ™2, «3 



Y^x^Yi^xip)) Yi 3 (x 3p ) 
L2 



n P ,9Zp) 



AM 



where 



7>(</>2p,</>3p) 



■ ,v n 2 ,n3 



P 



sin 



hp cos 



6 2 pSin £l+ ^+ 2 " 2 03pcos' 3 3 p x 



(cos20 2 p)-Prc3 2 (cos2</>3p) 



(10) 



where P^ ,fc are Jacobi polynomials and the coefficients AC„J'„ ' 1 3 are normalization factors. The quantity .fT = £1 + 
^2 + ^3 + 2(n 2 + n 3 ) is the so-called grand angular quantum number. The HH functions are the eigenfunctions 
of the hyperangular part of the kinetic energy operator. Another important property of the HH functions is that 
P K yf 1 ' 1 t2% l 2 n 2 n 3 (^p) are homogeneous polynomials of the particle coordinates of degree K. 

A set of antisymmetrical hyperangular-spin-isospin states of grand angular quantum number K, total orbital angular 
momentum A, total spin E, total isospin T, total angular momentum J, and parity 7r can be constructed as follows: 



KASTJtt 



12 

qKAZTJit 

p=\ 



(i,j,k,m) , 



(11) 
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where the sum is over the 12 even permutations p = k, m, and 



KASTJ71 



[s iSj ], 



(12) 



TT. 



Here, 3^"'^,'^f La „ 2 „ 3 (^p) is the HH state defined in @, and s, (t^) denotes the spin (isospin) function of particle 
i. The total orbital angular momentum A of the HH function is coupled to the total spin S to give a total angular 
momentum J and parity tt = (— 1)^1+^2+^3 _ The integer index fi labels the possible choices of hyperangular, spin and 
isospin quantum numbers, namely 



(j, = {£i,£ 2 ,£ 3 , L 2 ,n 2 ,n 3 , S a , S b ,T a ,T b } , 



(13) 



that are compatible with the given values of K, A, E, T, J and 7r. Another important classification of the states is 
to group them in "channels" : states belonging to the same channel have the same values of angular (ii, £ 2 , £3, L 2 ,A), 
spin (S a , Sb, E), and isospin (T a , T b , T) quantum numbers but different values of n 2 , n 3 . 

Each state \^kay,tj^ en t erm g the expansion of the 47V wavefunction must be antisymmetric under the exchange of 
any pair of particles. Consequently, it is necessary to consider states such that 



(14) 



which is true when the condition 



4 + S a + T a = odd 



(15) 



is satisfied. 

The number Mkahtj-k of antisymmetrical functions ■qt^ A ' STj7r having given values of K, A, E, T, J and tt but 
different combination of quantum numbers /i (see Eq. (^3j)) is in general very large. In addition to the degeneracy 
of the HH basis, the four spins (isospins) can be coupled in different ways to S (T). However, many of the states 
tyKAZTjir ^ p _ 1 _ _ ; Mkaztjtt are linearly dependent. In the expansion of a AN wave function it is necessary to 
include the subset of linearly independent states only, whose number is fortunately significantly smaller than the 
corresponding value of Mkayttj-k ■ 

The internal part of the wave function can be finally written as 



LSJJ z -k 
C 



E 

KAT.T 



E 



KAY.TJ-R 



(16) 



where the sum is restricted only to the linearly independent states. 

The main problem is the computation of the matrix elements of the Hamiltonian. This task is considerably simplified 
by using the following transformation 



4> 



KAT.TJ 



*(i,j;k;m) 



KASTJ >)$£ AST ^(1,2;3;4) . 



(17) 



The coefficients a^ J ^F Tj7T (p) have been obtained using the techniques described in Ref. OH} Then the kinetic energy 
operator matrix elements are readily obtained analytically, and the NN (3N) potential matrix elements can be 
obtained by one (three) dimensional integrals. The details are given in Ref. [3. 

The S'-matrix elements $ls l'S'(p) an0 - functions u M (p) occurring in the expansion of \^>^ SJJ ^ are determined by 
making the functional 



i^LS,L'S' 



°LS,L'S' 



(?) 



M 



N 



V6i 



us' JJ X 
1+3 



H — E3 — — 
2ji 



LSJJ 2 
1+3 



(18) 



stationary with respect to variations in the S~[ s L , s , and u^{p) (Kohn variational principle). Here E 3 is the 3 He ground- 
state energy. By applying this principle, a set of second order differential equations for the functions M^y^fp) is 
obtained. By replacing the derivatives with finite differences, a linear system is obtained which can be solved using 
the Lanczos algorithm. This procedure, which allows for the solution of a large number of equations, is very similar 
to that outlined in the Appendix of Ref. 03 and it will not be repeated here. 

The main difficulties of the application of the HH technique are the slow convergence of the basis with respect to 
the grand angular quantum number K, and the (still) large number of linearly independent HH states with a given K. 
Also a brute-force application of the method is not possible even with the most powerful computers available, so one 



6 



TABLE I: Quantum numbers of the first seventeen channels considered in the expansion of the state wave functions. See 
the text for details. 
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has to select a suitable subset of states [HS EH> ELi| • In the present work, the HH states are first divided into classes 
depending on the value of L = £\ + £2 + £3, total spin S, and 122, rc.3. In practice, HH states of low values of t\, £2, 
£3 are first included. Between them, those correlating only a particle pair are included first (i.e. those with ni = 0), 
then those correlating three particles are added and so on. The calculation begins by including in the expansion of the 
wave function the HH states of the first class CI having grand angular quantum number K < K\ and studying the 
convergence of a quantity of interest (for example, the phase-shifts) by increasing the value of K\. Once a satisfactory 
value of K\ = ifi max is reached, the states of the second class with K < K2 are added in the expansion, keeping all 
the states of the class CI with K\ = ifi, nax . Then K2 is increased until the desired convergence is achieved and so 
on. 

Let us consider, for example, the case J 71 ' = 0~ , where there is only one LS channel in the sum over L'S' of Eq. JSJ, 
namely L' = 1, S' = 1. Consequently, for this wave the S-matrix reduces to a value which is parametrized (as usual) 
as 6>£ S LS — ?yexp(2i(5), where 77 is known as the "elasticity parameter". Note that in the application of the Kohn 
principle given in (|18|l the value 77 = 1 is not guaranteed: it is achieved only when the corresponding internal part 
^jLSJJzTt j g we jj described by ^e HH basis. We have used the value of 77 as a check of the convergence of the HH 
expansion. In cases of poor convergence, the value of rj has been found to depend very much also on the choice of 
fL(yi), the function used to regularize the Coulomb function Gl- This function has been chosen to depend on a 
non-linear parameter 7, and thus a test of the convergence is performed by analyzing the dependence of rj on the 
parameter 7. At the beginning of the calculation, when the number of HH functions is not large enough to get 
convergence, rj is extremely dependent on the value of 7. The phase shift 5, however, depends less critically on 7. 
By increasing the number of HH components in the internal wave function we observe that 77 — ► 1, and that the 
dependence on 7 becomes negligible. Notice that the convergence rate has been found to depend on the value of 7; 
some critical values of this parameter exist where the convergence can be very slow. However, it is not difficult to 
find ranges of values of 7 where the convergence is fast and smooth and the final results are independent of 7. Since 
we are interested in the convergence of the HH function, we have chosen 7 in one of the "favorable" regions, where 
the convergence is achieved in a smooth and fast way. A detailed study on this subject will be reported elsewhere. 

Let us now briefly discuss the choice of the classes of HH states for the J T = CP case. Note that, since the wave 
under consideration is of negative parity, only HH functions with odd values of K and C have to be considered. 
Moreover, we consider in this work only states with total isospin T — 1, as the effect of the states with T — 2 should 
be negligible. The criteria used to select the appropriate classes of HH functions require that the states with lowest 
C be considered first. A few of the channels considered in the calculation have been reported in Table [I] The final 
choice of classes for the case J 71 " = 0~ is detailed below. 

1. Class CI. In this class are included the HH states with 712 = belonging to the channels 1 through 8 of Table [I] 
Note that the corresponding radial part of the HH functions depends essentially on cos</>3 P = r^j p and thus 
these states take into account two-body correlations (see Equation (TDJ)) . This part of the wavefunction is more 
difficult to construct due to the strong repulsion at short interparticle distances. 
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TABLE II: Summary of Measurements. The A y measurements were made using a proton beam and the cr{9) measurements 
were made using a He beam. The center-of-mass energies at which the two sets of measurements were made do not exactly 
agree due to differences in analyzing magnet calibrations for the two sets of experiments, and due to energy losses in the targets. 
All energies are in MeV. 



E c .m. 


A y measurements 


E p 




a (8) measurements 
E p 


E3 Ho 








0.74 


0.99 


2.97 


1.20 




1.60 


1.19 


1.59 


4.76 


1.69 




2.25 


1.67 


2.24 


6.69 


2.35 




3.13 


2.33 


3.11 


9.31 


3.03 




4.05 


3.02 


4.02 


12.06 



2. Class C2. This class includes HH functions belonging to the same eight channels as for class CI, but with 
n 2 > 0. Since cos</> 2p is proportional to the distance of particle k from the center of mass of the pair ij, these 
states, therefore, include also part of the three-body correlations. 

3. Class C3. This class includes the remaining T = 1 states of the channels having C — 1 (channels 9 through 17 
of Tabled). 

4. Class C4. This class includes the T = 1 states belonging to the remaining channels with £ = 3 and E = 1. 

5. Class C5. This class includes the T = 1 states belonging to the channels with C = 3 and E = 2. 

6. Class C6. This class includes the T = 1 states belonging to the channels with C = 3 and E = 0. 

7. Class C7. This class includes the T = 1 states belonging to the channels with L = 5. 

All states of the first four classes have a total spin E = 1 . The classification related to the total spin is important 
since we have observed that the component with E = 1 is the dominant one and requires more states to be well 
accounted for, while the E = 2 and E = components give only a tiny contribution to the phase shift (however, they 
are important for achieving 77 = 1). 

We have also calculated the phase-shifts of the states = + , 1 + , 2 + , 3 + , l - and 2~. The choice of the classes 
in these cases has been performed in the same way as discussed above. In order to test the accuracy reached by 
the theory, significant work has been done checking the convergence of the HH expansion in terms of the various 
classes. Some examples of the convergence for the phase shift 5 and elasticity parameter r\ using the AV18 potential 
are discussed in the Appendix. 



III. EXPERIMENTS 



Angular distributions of cross sections &{6) and proton analyzing powers A y were measured with high precision for 
proton— 3 He elastic scattering at several energies below 5 MeV. Silicon surface-barrier detectors, having an effective 
efficiency at these energies of 100.0 ± 0.1%, were used for both sets of measurements. In all of the measurements, a 
pulse generator signal was sent to the test input of the preamplifier of each detector and used as a measure of the 
dead-time of the data acquisition system. 



A. Cross-section measurements 



All of the measurements of the differential cross section were made using the TUNL supersonic gas-jet target, a 
refurbished and upgraded version of the target designed and built at the University of Erlangen-Niirnberg |55j . Target 
thicknesses of « 5 x 10 17 atoms/cm 2 were routinely obtained in a narrow jet (w 1 mm) for the measurements in this 
work, facilitating very reasonable counting times for the cross-sections measured. The target's lack of beam-degrading 
windows, well-defined geometry, and high-purity and high-density were ideal for these low-energy measurements. 

Since 3 He is a very expensive gas, the measurements were performed in inverse kinematics using a 3 He beam incident 
on a hydrogen gas-jet. Measurements of a(9) were made at the five 3 He energies listed in Table lll| Scattered 3 He 
particles can only be detected forward of #i a b ~ 19°. Detection of both the scattered 3 He and the recoiling proton 
allowed for the measurement of a(9) over a wide range of center of- mass- angles. 
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Beam 




FIG. 1: Scattering chamber setup for cr(8) measurements. The primary detectors are arranged symmetrically about the beam 
direction. For the relative u(0) measurements, the yields in the primary detectors were normalized to the counts in the fixed 
monitor detectors shown. The beam position and beam current were monitored using a set of slits behind the target. 

TABLE III: Collimation configurations for the relative cr(8) measurements. All lengths are in mm. All detectors were 16.6 cm 
from the target center, and fitted with a 5.1 cm snout. H is the horizontal collimator dimension and V is the vertical collimator 
dimension. The "Forward Angles" column shows the detector setups for forward angle measurements (with the monitors at 
55°) and the "Backward Angles" column shows the setup for backward angle measurements (with the monitors at 15°.) 







Forward Anj 


des 


Backward Angles 


He 


Detector 


Front 


Rear 


Front 


Rear 


(MeV) 




H x V (mm 2 ) 


H x V (mm 2 ) 


H x V (mm 2 ) 


H x V (mm 2 ) 


2.97 


Primary 


2.4 x 9.5 


0.8 x 9.5 


2.4 x 9.5 


0.8 x 9.5 




Monitor 


6.4 x 9.5 


6.4 x 9.5 


2.4 x 9.5 


0.8 x 9.5 


4.76 


Primary 


2.4 x 9.5 


0.8 x 9.5 


3.2 x 9.5 


1.6 x 9.5 




Monitor 


6.4 x 9.5 


6.4 x 9.5 


2.4 x 9.5 


0.8 x 9.5 


6.69 


Primary 


2.4 x 9.5 


0.8 x 9.5 


2.4 x 9.5 


0.8 x 9.5 




Monitor 


6.4 x 9.5 


6.4 x 9.5 


2.4 x 9.5 


0.8 x 9.5 


9.31 


Primary 


2.4 x 9.5 


0.8 x 9.5 


3.2 x 9.5 


1.6 x 9.5 




Monitor 


6.4 x 9.5 


3.2 x 9.5 


2.4 x 9.5 


0.8 x 9.5 


12.06 


Primary 


2.4 x 9.5 


0.8 x 9.5 


2.4 x 9.5 


0.8 x 9.5 




Monitor 


6.4 x 9.5 


6.4 x 9.5 


2.4 x 9.5 


0.8 x 9.5 



For each energy, a beam of 3 He ++ ions was accelerated with the FN tandem accelerator and then deflected onto 
the hydrogen gas-jet target. The elastically scattered protons and 3 He nuclei were detected by three pairs of primary 
detectors placed symmetrically about the beam direction, as shown in Fig. ^ Each detector was fitted with a pair 
of rectangular slits mounted in a cylindrical "snout." The dimensions of these slits are given in Table I1III The 
counts in each detector were normalized to the yield of scattered particles in a pair of monitor detectors also placed 
symmetrically about the scattering region. The angular range covered with the movable chamber detectors was 7.5° 
to 75°. When measuring forward angles the monitor detectors were placed at 55°, while for more backward angle 
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FIG. 2: (Color online) Typical energy spectrum of particles from the scattering of the 3 He beam on the hydrogen target. This 
spectrum was measured at #i a b = 15° at a 3 He beam energy of 6.73 MeV. The peak at channel 850 is caused by heavier gas 
contaminants. 



measurements the monitor detectors were placed at 15°. A cross-normalization was performed to maintain a consistent 
normalization for both sets of monitor detector angle settings. A typical spectrum is shown in Fig. [21 The beam 
position relative to the jet-target was monitored by a pair of horizontal and vertical slits behind the target region. 
The beam passing through the small slit opening was integrated by a Faraday cup which was electrically isolated from 
the slits. For each measurement the beam current was maximized on the Faraday cup, thus making sure that the 
scattering geometry remained the same from run to run. The currents from the slits and Faraday cup were summed 
and used to measure the total integrated charge on target. 

At lower beam energies, multiple scattering in the gas-jet decreased the measured total integrated charge and a 
small correction factor had to be applied. This factor was determined by frequently cycling the gas in the jet on and 
off and determining the effect of the gas presence on the integrated beam current. This was found to be a (6 ± 1)% 
correction at -E3 He = 2.97 MeV, reducing to a negligible correction at i?3 He = 12 MeV. 

The absolute normalizations of the u{9) measurements were performed with two methods. In the first method, 
proton- 3 He elastic scattering was normalized to 3 He- 40 Ar Rutherford scattering. Bombarding a gas-jet containing 
both hydrogen and argon with a 3 He beam, the ratio of the proton- 3 He yield to the 3 He- 40 Ar yield in the same 
detector at a given angle was measured. If the ratio Rt of hydrogen target thickness to argon target thickness is 
known, the ratio of p- 3 He counts to 3 He- 40 Ar counts at the same angle gives an absolute determination of a(9) at 
that angle. 

For measurements of the absolute normalization using this method, a small amount (~ 3%) of Argon was mixed 
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TABLE IV: Typical error budgets for both the mixed-jet and beam-switching absolute normalization methods. Statistical and 
systematic errors are listed. 





Mixed- Jet 




Source 


Type 


Error (%) 


Counting Statistics 


stat. 


0.3 


40 Ar/H 2 Ratio 


stat. 


0.6 




sys. 


0.5 


40 Ar - 3 He a{9) 


sys. 


0.8 


Angle Setting 


sys. 


1.0 




Beam- Switching 




Source 


Type 


Error (%) 


Counting Statistics 


stat. 


0.4 


Proton-proton a(E, 9) 


sys. 


0.7 


BCI Correction Factor 


sys. 


1.0 


Beam Energy 


sys. 


<0.1 


Angle Setting 


sys. 


<0.1 


Jet Reproducibility 


sys. 


0.8 


TABLE V: Overall systematic normalization error for each of the a(9) measurements. 


Equivalent E p [MeV] 




Error (%) 


0.99 




3.5 


1.59 




2.0 


2.25 




2.7 


3.11 




2.9 


4.02 




2.7 



with the hydrogen gas making the target jet. The gas was mixed by the manufacturer and the Ar/H2 ratio Rt 
was determined to an accuracy of 2% by gas chromatography 56]. Rt was also measured using a proton beam at 
E p = 2.24 MeV, at angles where proton- 40 Ar elastic scattering is known to be described by the Rutherford formula 
within a percent, as calculated using several different sets of optical model parameters |57|. Using the well-known 
proton-proton elastic-scattering cross-section, which was obtained from the high-accuracy phase-shift analysis of the 
Nijmegen group |58l l59|. determinations of Rt using this method agreed within error with the gas-chromatography 
measurements . 

Similarly, this mixed-jet method also relies on 3 He- 40 Ar elastic scattering being described correctly by the Rutherford 
scattering formula. Using the DWBA code DWUCK4 and optical model parameters from Ref. |6lJ, it was determined 
that the 3 He- 40 Ar elastic scattering cross-section is within 5% of the Rutherford prediction out to #i a b = 40° for the 
three lowest 3 He energies in Table ITT1 

The other technique for determining the absolute normalization of the cr{9) angular distributions was a beam- 
switching method, in which the product of detector solid-angle and target-thickness was determined using a proton 
beam incident on a hydrogen jet, via the known proton-proton elastic scattering cross section. A 3 He beam at the 
proper energy was first tuned onto the hydrogen jet target; this was followed by irradiating the jet with a proton 
beam with the same magnetic rigidity. This procedure allowed for minimal beam-transport adjustments; the source 
inflection magnet before the tandem accelerator and the accelerator terminal potential were adjusted so that both 
beams passed into the chamber with the same beam tune. This ensured the beam-target geometry was the same for 
each beam. Scattered particles from each beam were detected by three pairs of fixed-angle detectors. This procedure 
was repeated several times to ensure reproducibility. 

The beam-switching technique was used at energies at which the mixed-jet method was not feasible. Both techniques 
were used at several energies, as a cross-check, and the results from both methods agreed within errors. Typical error 
budgets for each method are shown in Table iTVl and the overall systematic normalization errors are listed in Table 

m 
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FIG. 3: (Color online) Chamber setup for A y measurements. 



TABLE VI: Horizontal (H) and vertical (V) detector slit dimensions for A y measurements. The front and back collimators for 
each detector were the same dimensions, as listed. The distance from the target center to the detector face R, and the length 
of the detector snout S are also listed. The detector pairs were spaced at 15° intervals. 







Forward Pair 






Middle Pair 




Backward Pair 




E p 


H x V 


R 


S 


H x V 


R 


8 


H x V 


R 


S 


(MeV) 


(mm 2 ) 


(cm) 


(cm) 


(mm 2 ) 


(cm) 


(cm) 


(mm 2 ) 


(cm) 


(cm) 


1.60 


1.6 x 9.5 


12 


6.4 


1.6 x 9.5 


12 


5.1 


2.4 x 9.5 


14 


5.1 


2.25 


1.6 x 9.5 


10.2 


5.1 


1.6 x 9.5 


12 


5.1 


2.4 x 9.5 


14 


5.1 


3.13 


1.6 x 9.5 


10.2 


5.1 


1.6 x 9.5 


12 


5.1 


2.4 x 9.5 


14 


5.1 


4.05 


1.6 x 9.5 


10.2 


5.1 


1.6 x 9.5 


12 


5.1 


2.4 x 9.5 


14 


5.1 



B. Analyzing power measurements 

The measurements of A ,, w ere made utilizing the atomic beam polarized ion source at TUNL |62j | via a two 
polarization state method |63j with fast state switching 64] . This polarized proton beam was accelerated to the 
desired energy with the tandem accelerator and tuned onto a gas-cell target inside the 62 cm diameter scattering 
chamber. The gas-cell, employing a 2.3 /im Havar foil, was filled with 1 atm of 3 He gas, and was mounted on the 
target-rod which was supported from the top of the chamber. This allowed the cell to be raised, allowing the beam 
to directly enter the polarimeter. A schematic of the experimental setup for the A y measurements is shown in Fig. [3J 
and the collimation setup is detailed in Table IVll These collimators limited the view of the detectors to only protons 
scattered from 3 He gas in the cell and not those scattered from the cell entrance and exit foils. A y data was taken 
only at angles for which foil-scattering was negligible. 

Since there is significant energy loss (particularly at the lowest energies) in the cell foil, the incident beam energy 
was adjusted so that the beam reached the desired energy at the cell center. The energy losses in the cell foil and 3 He 
gas were modeled by the computer code SRIM-2000 J65|. The proton energies at the center of the gas-cell are listed 
in Table HD 

The polarization of the proton beam was monitored on-line with a polarimeter based on 4 He(p,p) 4 He clastic 
scattering jggj. Periodically during the experimental runs the beam energy was raised either to 6 MeV or 8 MeV, 
where the analyzing power for the polarimeter is very close to unity. This was done once every two to three hours. 
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The polarization state of the beam was switched approximately three times a second and was typically « (67 ± 2)%. 
A 2% systematic error on the A y measurements arises from the error in beam polarization determinations. 

IV. COMPARISONS WITH THEORY 

In this section, the experimental data are presented and compared with the results of the theory reviewed in Sec, till 
The results for the differential cross sections and analyzing powers are presented in Section TlV Al and Section TlV Bl 
respectively. Note that the a (9) data are designated by their equivalent proton lab energy E p , despite the data 
being taken in inverse kinematics. Finally, the theoretical predictions of other observables reported in Ref. 1671 are 
presented and discussed in Section llV CI The calculations presented were performed using the Argonne vis [47| NN 
potential (AV18 model), and with the vm NN potential with the inclusion of the Urbana IX 3N force (AV18/UIX 
model) . The corresponding phase shift and mixing angle parameters calculated with the HH expansion have reached 
a noticeable degree of convergence, as discussed in great detail in the Appendix. 

A. Differential Cross-Section cr(0) 

The measured differential cross sections uiff) at the five energies considered here are presented and compared with 
the existing data j2^| in Fig. 0] The results of the previously described calculations for the AV18 potential (dashed 
lines) and AV18/UIX model interactions (solid lines) are also shown. 

When comparing the data of this work with previous measurements of a{9) at the same energies, the agreement 
(in general) is quite good; there is a marked improvement in the size of the error bars, and the new data sets contain 
many more data points. The precision is much better than the data of Reference [2^, slightly better than that of 
Reference |68|. and is comparable to that of Reference j2^. At E p = 4.02 MeV, there is good agreement between 
the current data and those of Reference 68] but with slightly smaller error bars. There is no previous data known to 
exist at E p = 3.11 MeV. 

As can be seen in Fig. there is a general agreement between the theoretical and the experimental results. 
For small angles, the cross section is dominated by the Coulomb scattering and a good agreement is observed (the 
exception being two points at small 9 c . m , at E p — 4.02 MeV). At 9 c . m , — 90°, the contribution of the L = 1 waves 
vanishes, and therefore the cross section is almost completely due to the L = phase-shifts. As discussed in the 
Appendix, there are no problems in the calculation of the L = phase-shifts from the numerical point of view, and 
therefore cr(90°) is an unambiguous test of the underlying nuclear dynamics. We observe a sizeable 3N force effect in 
this region (the minimum), which tends to decrease as E p is increased. This is consistent with the increased binding 
of the 3 He when the 3N force is included. 

As 9 c . m , approaches 180°, the predicted cross section becomes quite sensitive to the L = 1 phase-shifts (states with 
J n = - , 1~, and 2~). The calculations slightly underestimate the cross sections there, particularly when the 3A" 
force is included. This problem is somewhat analogous to that found in n- 3 H elastic scattering in the peak region 
at E n w 3 MeV [l4|, as mentioned in Section ITAl In that case, the calculations based on the standard NN and 3N 
forces (as used here) are found to underestimate the total cross section by 20% on the peak [l4|. These problems 
probably arise from an incomplete knowledge of either the NN or the 3N interaction in P-waves, and therefore are 
closely related to the N-d A y puzzle. This becomes more evident in the study of the p- 3 He analyzing powers presented 
below. 



B. Proton Analyzing Power A y 

The measured proton analyzing power A y at the four energies considered here are presented and compared with the 
existing data in Fig. EJ Additionally, at E p = 1.0 MeV the experimental data of Ref. [22 are shown. The calculations 
obtained with the AV18 (dashed lines) and AV18/UIX (solid lines) models are also shown. Note that A y steadily 
grows as EL is increased. There is good agreement between the new measurements and the older ones reported in 
Refs. [1(JlI3- Note, however, that the present measurements are noticeably more precise, in particular at E p = 1.60 
and 2.25 MeV. 

The theoretical calculations clearly underestimate the data at all energies. The 37V force of the Urbana-type 
has a very little effect at low energies, but its effects are larger at E p = 3.13 and 4.05 MeV. They are, however, 
clearly insufficient to resolve the discrepancies with the data. The present results confirm the disagreement previously 
reported in Refs. [Tolll2| . A plot of the relative difference between experiment and theory at the maximum A y value in 
the angular distribution as a function of proton energy is shown in Fig. [|J] Note that this difference is nearly constant 
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FIG. 4: (Color online) The measured p- 3 He elastic differential cross sections (solid circles) at five different energies are compared 
with the data reported in Ref. 23> (open squares), Ref. [22 (open diamonds), and Ref. [63 (open circles). The curves show the 
results of the theoretical calculations for the AV18 (dashed lines) and AV18/UIX (solid lines) potential models. 



as the energy is changed. This is similar to what is observed in N-d scattering, though the discrepancy in the p- 3 He 
case is about 50% larger. The A y observable is very sensitive to the L = 1 phase-shifts, and in particular to the 
combination of phase shifts A = S( 3 P 2 ) - [S( 3 Pi) + <5( 3 P )]/2 [hJ. The value of A is predicted (using the AV18 and 
AV18/UIX models) to be smaller than the one extracted from the data. It is interesting to note that this is analagous 
to the N-d case, in which the splitting between the 4 Pi/2 phase and the average of the 4 -p3/2 and 4 P$/2 phases is 
too small to reproduce the observed A y . It would be very interesting to see if new terms in the 37V interaction could 
explain both the N-d and p- 3 He A y discrepancies. Work in this direction is in progress. 
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FIG. 5: (Color online) The measured p- 3 He proton analyzing power A y (solid circles) at five different energies are compared 
with the data of Ref. (open squares), Ref. [52 (open diamonds), and Ref. [63 (open circles). The curves show the results of 
the theoretical calculations for the AV18 (dashed lines) and AV18/UIX (solid lines) potential models. 



C. Other Observables at E v = 4.05 MeV 

At E p — 4.05 MeV, measurements of other p- 3 Hc observables (the spin correlation coefficient A yy and the 3 He 
analyzing power Aq v ) exist [6?| . The comparison between the results of the present calculation and these data are 
shown in Fig. The measurements have rather large error bars, and no clear conclusion about the agreement with 
theory for the A yy measurements can be reached. However, Aq v does appear to be under-predicted at the maximum. 
Indeed, Aq v is particularly sensitive to the combination <5( 3 i"2) — <5( 3 -Pi)- On the other hand, the observable A yy is 
quite sensitive to e(l + ), the mixing parameter of the J T = 1 + state. More precise measurements of these observables 
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FIG. 6: (Color online) Relative difference as a function of energy between the theoretical predictions and measured values for 
A y at the peak of the distribution for p- 3 He scattering (•). Also shown are results from Ref. j3f° r P~d C) an d n-d scattering 
(❖)■ 

could provide much-needed input for an experimental phase-shift analysis and therefore permit a better understanding 
of the discrepancy with theoretical calculations. Such measurements are currently underway [69^. 

V. CONCLUSIONS 

In this paper the solution of the Schrodinger equation for four-nucleon scattering states has been obtained using the 
HH function expansion. The main difficulty when using the HH basis is its large degeneracy. Accordingly, a judicious 
selection of the HH functions giving the most important contributions has been performed. For this work, the HH 
functions have been divided into classes, depending on the number of correlated particles, the values of the orbital 
angular momenta, the total spin quantum number, etc. For each class, the expansion has been truncated so as to 
obtain the required accuracy. The study of the convergence of p- 3 He elastic scattering phase-shifts and observables 
reported in the Appendix has shown that good accuracies are achieveable and a powerful method to extrapolate the 
results has been also discussed. When applied previously for n- 3 H elastic scattering, the HH method has been proved 
to give results in good agreement with other theoretical techniques [T^| . 

We also reported measurements of the cross-section a{&) and proton analyzing power A y for p- 3 He elastic scattering 
over the range of energies 1.6 MeV < E p < 4.05 MeV. Additionally, a (8) measurements were obtained at E p = 0.99 
MeV. Analyzing powers were large and increased in magnitude by more than a factor of 3 as the energy was increased 
from 1.6 to 4.05 MeV. Both the A y and a (6) measurements have higher statistical precision at more angles, and 
smaller and well-understood systematic errors than those existing previously. 

There is good agreement between the cross section data and the calculations when the 3N potential is included. 
However, there are large differences between the A y data and theory. At the maxima of the A y angular distribution, 
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FIG. 7: (Color online) 3 He analyzing power Ao y and spin-correlation coefficient A yy at E p — 4.05 MeV. Data are from Ref. loll 
The curves show the theoretical predictions from the AV18 (dashed lines) and AV18/UIX (solid lines) potential models. 



the theory underpredicts experimental values by about 50 %. The inclusion of the 3iV potential produces only a small 
percentage change in the predicted analyzing powers and hence has little effect on the magnitude of the disagreement. 
This disagreement is remarkably similar to (and twice as large as) the " A y puzzle" observed for nearly 30 years for 
N-d scattering. 

The present calculations were extended to include A yy and Ao y for which there are measurements at 4.05 MeV fr?^ . 
The inclusion of the 3iV force has some influence on predictions of both Ao y and A yy . The calculations for these two 
observables are much closer to the experimental data although the data have much larger errors than for A y . More 
precise measurements of Ag y and A yy could help to define the phase shifts a nd p rovide a better understanding of the 
origin of this new "A y puzzle." Such measurements are currently underway 69] . 



APPENDIX: CONVERGENCE OF THE CALCULATED PHASE-SHIFTS 



In this Appendix, we discuss the convergence of the calculated phase-shifts. At the energies of interest here, p- 3 He 
scattering is dominated by the L = and 1 waves. The convergence of the HH expansion of -qi^ SJJ '' n: for the L = 
waves (J w = + , 1 + ) can be obtained by including a rather small number of channels. This is due mainly to the 
Pauli principle which limits overlaps between the four nucleons. As a consequence, the internal part is rather small 
and does not require large number of channels to be well described. 

On the other hand, for L = 1 waves (J" = Q~, 1~ and 2~) the convergence rate is slow and many channels have to 
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TABLE VII: Convergence of inelasticity parameter r\ and phase-shift 5 at E v = 


4.05 MeV corresponding 


; to the inclusion 


in the internal part of the 


wave function of the different classes CI 


- C7 in which the HH basis has been subdivided. The AV18 


potential is 


considered here with the inclusion of the point Coulomb interaction. 






Kr 


K 2 


Ks K 4 K 5 




V 


8 (deg) 


21 








1.00032 


10.649 


31 








1.00069 


11.484 


41 








1.00107 


11.882 


51 








1.00133 


12.060 


61 








1.00146 


12.136 


61 


11 






1.00139 


12.599 


61 


21 






1.00131 


12.897 


61 


31 






1.00134 


13.020 


61 


37 






1.00136 


13.055 


61 


37 


11 




1.00064 


15.284 


61 


37 


21 




1.00049 


15.923 


61 


37 


31 




1.00048 


16.105 


61 


37 


35 




1.00048 


16.132 


61 


37 


35 11 




1.00045 


16.256 


61 


37 


35 21 




1.00040 


16.646 


61 


3 1 


35 2o 




1.00040 


16.727 


61 


37 


35 31 




1.00040 


16.794 


61 


37 


35 31 3 




1.00012 


16.877 


61 


37 


35 31 7 




1.00002 


17.003 


61 


37 


35 31 11 




1.00000 


17.101 


61 


37 


35 31 15 




1.00000 


17.157 


61 


37 


35 31 19 




1.00000 


17.191 


61 


37 


35 31 19 


11 


1.00000 


17.194 


61 


37 


35 31 19 


11 11 


1.00000 


17.219 



be included. In these cases, the interaction between the p and 3 He clusters is very attractive (it has been speculated 
that some AN resonant states exist) and the construction of the internal wave function is more delicate. 

Finally, the contribution from L = 2 waves is rather tiny, since the centrifugal barrier does not allow the two 
clusters to come close, and the corresponding phase-shifts can be calculated with good approximation by neglecting 
the internal part i^^ SJJ ^ _ Contributions from L = 3 or higher waves has been disregarded since they are assumed 
to be negligible. 

Let us discuss in detail the convergence of the HH calculation of the CP phase-shift; the other J 77 states will be 
reported elsewhere. As shown in Section^ for this state L, S = 1, 1 and n can be parametrized as 77exp(2i<5). 
The results obtained for 77 and 6 at E p = 4.05 are reported in Table IVTT1 Here we have considered the AV18 potential 
model 47); however, the electromagnetic interaction has been limited to just the point-Coulomb potential. We have 
used 1/Mjv = 41.47108 MeV fm 2 . We study the convergence as explained in Section[Hj and the results presented in 
Table IVHI are arranged accordingly. For example, the phase-shift S reported in a row with a given set of values of 
K\ , . . . , Kj has been obtained by including in the expansion all the HH functions of class Ci with K < Ki, i = 1, . . . , 7. 

The convergence of the class CI is rather slow and a fairly large value of K has to be used. The inclusion of the 
second and third classes increases the phase-shift by about 4°. The class C4 contributes for additional 0.6°. The 
number of the states of this class increases very rapidly with K4 but fortunately the convergence is reached around 
K = 21. Up to now, the expansion includes only states with S = 1. The contribution of the states with E = 2, first 
appearing when the class C5 is considered, is rather small, and the contribution of the states with E = (class C6) 
is practically negligible. The contribution of the class C7 is also small. Since the number of states of this class is 
very large (there are 121 channels with l\ + £2 + (3 — 5 for J 71 ' = _ ) when confronted with a very tiny change of 
the phase-shift, a selection of the states has to be performed to save computing time and to avoid loss of numerical 
precision. In the present example, only the channels with (^1,^2,^3) = (1, 2, 2) have been found important. Note that 
at lower energies the convergence is noticeably faster (see below) . 

The convergence rate when considering the AV18/UIX interaction model is similar to the AV18 case. Since the 
models most frequently used for the 3iV interactions lack a strongly-repulsive core at short interparticle distances, the 
convergence rate of the various classes is found not to change appreciably. 

In order to obtain a quantitative estimate of the "missing" phase-shift due to the truncation of the HH expansion 
of the various classes, let us consider 8{K\, K 2 , K 3 , K4, K 5 , K 6 , if 7), the phase-shift obtained by including in the 
expansion all the HH states of the class CI with K < K\, all the HH states of the class C2 having K < K2, etc. Let 
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TABLE VIII: Convergence of CP phase-shift at Ep = 1.00, 2.25 and 4.05 MeV corresponding to the inclusion in the internal 
part of the wave function of the different subsets of HH basis. The AV18 potential is considered here with the inclusion of the 
point Coulomb interaction. In the last row, the total missing phase-shifts computed as described in text have been reported 
(we have estimated e~ 2 ~' w 0.8). 



E p = 1.00 MeV E p = 2.25 MeV E p = 4.05 MeV 



Ki 


K 2 


K 3 


Ki 


K 5 


K 6 


K 7 


S (deg) 


6 (deg) 


S (deg) 


51 


27 


25 


21 


1 


1 


1 


1.867 


7.401 


16.364 


53 


29 


27 


23 


3 


3 


3 


1.880 


7.467 


16.559 


55 


31 


29 


25 


5 


5 


5 


1.891 


7.526 


16.724 


57 


33 


31 


27 


7 


7 


7 


1.902 


7.583 


16.878 


59 


35 


33 


29 


9 


9 


9 


1.912 


7.632 


17.010 


61 


37 


35 


31 


11 


11 


11 


1.919 


7.668 


17.106 








5« 








0.028 


0.144 


0.384 



us compute 

Ax(K) = 8(K,0,0,0,0,0,0)-S(K -2,0,0,0,0,0,0) , (A.l) 

A 2 (K) = S(K\,K, 0,0, 0,0,0) - 5(K\, K- 2, 0,0, 0,0,0) , K\ = 61 , (A.2) 

A 3 {K) = 5{K[,K^, K, 0,0, 0,0) - S(R\,K^, K -2, 0,0, 0,0) , K~ 2 = 37 , (A.3) 

and so on. The values obtained for Aj, i = 1, 5 are shown in Fig. [SJfor the E p = 4.05 MeV on a logarithmic scale. 
As can be seen, all the differences Ai through A5 decrease exponentially, and approximately with the same decay 
constant. For a given K, however, there is a clear hierarchy Ai(K) ^> A 2 (K) ss A 3 (K) « A4(K) ^> As(K). Note 
that there are slight fluctuations in A(K) as K is increased (this is evident particularly for A5). The phase-shift 
differences for the classes C6 and C7 are not shown since they are tiny. 

From the simple behavior observed in Fig. |SJ we can readily estimate the missing phase-shift due to the truncation 
of the expansion to finite values of K = K. Let us suppose that the states of class i up to K = K have been included 
and used to compute Ai(K). From the observed behavior Ai(K) oc exp(— 7 AT), the "missing" phase-shift 8^ due to 
the states with K = K + 2, K + 4, can be estimated as 

6? = c( 7 ) MK) , (A.4) 

where 

00 

c( 7 ) - y e -7(*-^ = JL- , and x = . 

t * y — % 

K=K+2,K+4,... 

For example, consider the missing phase shift for the class CI. For K = 61, A\(K) = 0.009° and x ~ 0.8. Therefore, 
Sf 1 ps 0.04°, a rather small quantity. The missing phase-shift of the other classes can be estimated in the same way. 

However, to estimate the total missing phase-shift 5^ due to the truncation of the expansion of the first class up to 
K < Ki, of the second class up to K < K 2 , etc., we cannot simply add the A^ 1 , i = 1, . . . , 7 so obtained. The inclusion 
of the HH states of classes C2, C3, . . ., also alters the convergence of class CI by a small amount, etc. To study the 
"full" rate of convergence, we have taken advantage of the fact that the various A(K) show a similar convergence 
behavior (with approximately the same decay constant 7) and that the coefficient c(j) denned in Eq. (|A.4(i does not 
depend on K. Let us then consider 

A T {K X , ...,K 7 ) = &{K U ...,K 7 )- 5{K X - 2, . . . , K 7 - 2) . (A.5) 

From the above discussion, we can estimate the total missing phase shift as 

S¥ = T^— A T (K U ...,K 7 ) , x = e-^. (A.6) 

1 — x 

As an example, the values for Ai{K\, . . . , K 7 ) computed at E p = 1.00, 2.25, and 4.05 MeV are reported in Table IVlTTl 
from which it is possible to derive the values of 7 and then of S^ 1 . The computed values of 6jf using Eq. I A. 61 are 
reported at the bottom. 

As can be seen in Table IVHll 5^ is estimated to be rather small at E p — 1.00 and 2.25 MeV. However, for the 
largest energy the convergence seems not to be completely under control and higher values of K\ through K 7 should 
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FIG. 8: 0~ phase-shift differences for p- 3 He elastic scattering at E p — 4.05 MeV for the classes CI (circles), C2 (squares), C3 
(up triangles) C4 (asterisks), and C5 (down triangles) as function of the grand angular value K (see the text for more details). 
The potential used is AV18. 



be employed. In any case, we can see that the missing phase-shift is less than 2 %. Analogous problems have been 
found for the 1~ and 2~ states, whereas for the other states the convergence did not present any difficulty. 

Since at the moment the inclusion of a greater number of states would require a significant increase in computing 
time, we have preferred to use the extrapolation outlined above for obtaining estimates for the converged phase-shifts 
and mixing parameters for the J 71 " = CP, l - and 2~ states. To show the effect of the extrapolation on the observables, 
we present in Fig.[5]the results for four p- 3 He observables at E p — 4.05 MeV and calculations using the AV18 potential. 
The dashed and thin solid curves have been obtained using the 0~ phase shift calculated with different values for 
K\, . . . ,Ky. More precisely, the dashed (solid) curves have been obtained using the value S = 16.364° (17.106°) 
obtained with the choice of Ki, . . . , Kj reported in the first (sixth) row of Table IVTTT1 The thick solid curve has been 
obtained using the extrapolated value 5(61, 37, 35, 31, 11, 11, 11) + « 17.5°. The other phase-shift were taken with 
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FIG. 9: (Color online) Four p- 3 He elastic scattering observables at E p = 4.05 MeV calculated using different values for the 
CP phase-shift. The dashed and solid curves were obtained using different truncations of the HH expansion corresponding to 
choices of (Ki, . . . , Kr) reported in the first and sixth rows of Table IVTTT1 The thick solid curves have been obtained using 
the extrapolated value for this phase-shift, obtained as explained in the text. The experimental data are from Ref. Ir3^| (open 
circles) and from the present work (solid circles). The potential used is AV18. 



their final values (in particular, for the 1~ and 2~ phase-shifts, we used the extrapolated values obtained using a 
similar procedure as described above). The four observables considered in the Figure are: the differential cross section 
<j{9), the proton analyzing power A y , the 3 He analyzing power Ag y , and the spin correlation coefficient A yy . 

As can be seen, there is a good convergence for a (6), A yi and A yy . The observable A 0y is more sensitive to S, and 
the convergence is more critical. In any case, however, the thick solid curves are rather close to the thin solid curves, 
showing that the convergence has been nearly reached. A similar analysis has been performed also for the 1~ and 2~ 
phase-shifts and mixing parameters, with similar findings. Therefore, we can conclude that the convergence of the 
HH expansion is sufficiently good to obtain nearly correct predictions for the p- 3 He observables and allowing therefore 
for meaningful comparisons between theory and experimental data. The calculated phase shifts arc in reasonable 
agreement with those calculated by other methods (for example, those in Ref. IT3) . 
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